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Abstract 

Using the quark eigenmodes computed on the lattice with the overlap-Dirac operator, we inves¬ 
tigate the spatial distribution of the chiral condensate around static color sources corresponding 
to quark-antiquark and three-quark systems. A flux structure of chromo fields appears in the 
presence of such color charges. The magnitude of the chiral condensate is reduced inside the color 
flux, which implies partial restoration of chiral symmetry inside hadrons. Taking a static baryon 
source in a periodic box as a toy model of nuclear matter, we estimate the magnitude of the chiral 
symmetry restoration as a function of baryon matter density. 
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I. INTRODUCTION 


The low-energy dynamics of Quantum Chromodynamics (QCD) is characterized by two 
prominent properties, i.e. chiral symmetry breaking and conhnement. In the QCD vacuum, 
chiral symmetry of quark fields is spontaneously broken, as probed by an order parameter 
(gg), the vacuum expectation value of the scalar density operator gg. Through the Banks- 
Casher relation |1], this chiral condensate S = — (gg) can be related to the spectral density 
of the low-lying Dirac eigenvalues. In the presence of valence quarks, or color sources, some 
modihcation of the vacuum is expected around them. Strictly speaking, it is no longer the 
vacuum, i.e. the lowest energy state of the system, but we use this terminology having in 
mind an application to the study of finite density QCD. Although the static color sources 
as considered in this work, is only a crude approximation to the hnite density system, the 
study of the vacuum modification may shed light on the states of finite density QCD, which 
has been a subject of active research (see, for instance m)- 

The other interesting property of low-energy QCD is the confinement of quarks, which 
is characterized by the linearly rising potential between static color sources. Putting a pair 
of static quark and anti-quark in the vacuum, a color flux-tube emerges between them and 
leads to a linear increase of the energy as a function of the separation. This flux-tube 
structure has been observed in lattice QCD calculations by monitoring the action density or 
chromo-electric (or chromo-magnetic) field |3HS]. We expect that such flux-tube structure 
is reflected in the low-lying fermion eigenmodes, because the Dirac eigenmodes carry the 
information of their background gauge field conhguration. Indeed, the QCD held-strength 
tensor can be reconstructed using the fermion eigenmodes [6]. 

In this paper we present a lattice study of the spatial distribution of the chiral condensate 
in the presence of static color charges. We consider quark-antiquark and three-quark systems 
represented by Wilson loops to mimic the mesonic and baryonic states, respectively. We use 
the lattice data of the Dirac eigenmodes calculated on the gauge configurations generated 
with 2-|-l-flavor dynamical overlap fermions [7]. With the overlap fermion formulation [8l|9], 
chiral symmetry is exactly realized on the lattice, which is important in the study of the 
low-lying Dirac eigenmodes, as they are very sensitive to any small violations of chiral 
symmetry. The lattice data used in this work have this nice property, and indeed were 
successfully applied to the extraction of the chiral condensate in the vacuum [TUHIS]- 
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The organization of this paper is as follows. In Section |TT| we describe the method to 
constrnct the local chiral condensate qq{x) by nsing the overlap-Dirac eigenmodes, and show 
its distribntion in the vacnnm. In Sections |III| and [TV] we investigate the spatial distribntion 
of the local chiral condensate aronnd the static color sonrces. Section El is devoted to a 
snmmary. Preliminary reports of this work are fonnd in [miTb] . 

II. TOPOLOGICAL STRUCTURE OF THE QCD VACUUM 

We investigate the topological strnctnre of the non-pertnrbative QCD vacnnm in terms 
of the eigenmodes of the overlap-Dirac operator. It preserves exact chiral symmetry, and 
the relation to the topological charge of backgronnd gange field confignration is manifest, 
i.e. the index theorem, at least for smooth enongh backgronnds im. 

In this paper, we nse the 2-|-1-flavor dynamical overlap-fermion confignrations gener¬ 
ated by the JLQCD Collaboration [7j. Their lattice volnmes are 16^ x 48 and 24^ x 48 
at a single inverse lattice spacing a~^ = 1.759(10) GeV. The dynamical qnark masses are 
mud = 0.015a“^ and = 0.080a“^. The global topological charge is fixed at Q = 0 to 
avoid the problem of divergent molecnlar dynamics force in the simnlations |T8] . It indnces 
hnite volnme effects [12], which would not be significant for the relatively local observables 
considered in this study. Most of the results are obtained on the larger lattice (24^ x 48) 
where the number of independent conhgurations is 50. 

In the following we describe the profile of the low-lying eigenmodes on these lattices. 

A. Local chiral condensate qq{x) 

The massless overlap-Dirac operator is given by nu 

Dov(O) = mo [1 + 75Sgn Hwi-mo )], (1) 

with the hermitian Wilson-Dirac operator Hw{—mo) = Here, sgn denotes the 

matrix sign function. Introducing the quark mass m^, the overlap-Dirac operator is modihed 
as 

^ov(mq) = (^1 “ ^ov(O) mg. (2) 

This form cancels 0{a) discretization effects, together with a proper rotation of the fermion 
helds in the observables. 
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We define the eigenfunction 'il)\{x) associated with an eigenvalue A of the massless overlap- 
Dirac operator 

= X'lpxi.x), (3) 

where the eigenfunction 'ip\{x) is normalized as ~ 1- Using them we may 

expand the “local chiral condensate” qq{x) in terms of the eigenmodes, i.e. 

i^l{x)i^x{x) 


qq(x) = - ^ 


(4) 


for a valence quark mass rUg. This relation represents a self-contracting fermion loop contri¬ 
bution from and to the scalar density operator. If the measured observables do not include 
other light quark helds to be contracted, then the substitution Q is justihed. The correla¬ 
tion functions of qq{x) with the Wilson-loop are in this class of observables. 

The chiral condensate (qq) is given by an ensemble average of qq{x) without insertions 
of other operators. By averaging over space-time, this quantity is simply written in terms 
of only the eigenvalues because of the normalization condition for 'ipx{x). Thus the relation 
between the chiral condensate and the spectral density p(A) of the Dirac eigenvalues is 
established. In the chiral limit it reads S = 7rp(0), i.e. the Banks-Casher relation [I]. 


B. Action and topological charge densities in terms of the Dirac eigenmode 

Since the gauge held strength tensor is dehned through the covariant derivative 
as = [D^, D^], it can also be related to the Dirac operator [B]. Here we briehy reproduce 
the derivation. 

The square of the Dirac operator Ip = is decomposed as 

[pi.x)f = J2 dI{x) + (5) 

By multiplying 7^7^ and taking a trace with respect to the Dirac indices, the held strength 
tensor is expressed as 

Ff,u{x) = -^tr . (6) 

Therefore, by expanding the Dirac operator in terms of its eigenvectors 'pxix), an expansion 
of the held strength is obtained: 

Ff,u{x) = '^X‘^f^u{x)x, ft,u{x)x = ^'Pl{xht,-fui^x{x). (7) 

A 
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FIG. 1: Eigenvalue (in MeV) of the low-lying eigenmodes on the 24^ x 48 lattice. By including 
240 eigenmodes, we can cover the range of A < 300 MeV. 


Using this decomposition, the action and topological charge densities are expressed as 

p{x) (8) 

A,A' 

qtop{.x) = tT^'^X^X'‘^f^u{x)\ff,u{x)x', (9) 

A,A' 

respectively. Here, trc denotes the trace with respect to the color indices and f^i,{x)x = 

2^IJ.i'pu fpfri.x') X- 

So far the expressions are exact, but in the numerical studies we introduce a truncation of 
the summation over the eigenmodes. This truncation acts as a hlter to cut UV fluctuations 
above Amax- On the ensembles of 16 ^ X 48 and 24^ X 48 lattices, we calculated 160 and 
240 pairs of eigenvalues and eigenvectors of Dov, respectively. Then, the eigenvalues after 
correcting the 0(a) effect, Im A/(l — Re A/2mo), cover the region between ±300 MeV, as 
shown in Figure In the measurements of the correlation between qq{x) and the Wilson 
loops, we monitor the dependence on the number N of the eigenmodes included and conhrm 
that the results saturate at least above 200 MeV. Some examples will be shown later. 

Before presenting the results, we show some snapshots of the eigenmodes. The index the¬ 
orem dictates that exact zero-modes are associated with topological excitations of the gauge 
held. This suggests that the near-zero modes are superpositions of such local topological 
objects. Using a truncation at iV = 20, we visualize the low-mode contributions to the local 
chiral condensate qq{x) (|^, the action density p{x) (|^ and the topological charge density 
qtop{x) ([^ in the panels (a), (b) and (c) of Figure]^ respectively. They show tomographic 
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FIG. 2: Snapshots of (a) the local chiral condensate, (b) the action density and (c) the topological 
charge distributions observed with the sum of 20 lowest-lying eigenmodes. These pictures show the 
same T-X slice of a 24^ x 48 lattice on a representative gauge configuration. The chiral condensate 
local fluctuations are correlated with the local topological charge measurements. 


images on a certain T-X slice of the four-dimensional lattice extracted from a given gauge 
configuration of size 24^ x 48. 

As one can see in Fig. (a), the local condensate qq{x) forms a cluster structure. At the 
same space-time points of the cluster, the action density shows peaks (the panel (b)). More 
importantly, the topological charge density has positive and negative islands stretching over 
several lattice spacings at the same space-time points. Such observation is not new; indeed, 
there are lattice studies using the overlap-Dirac operator [201 El] showing the similar prohle 
of the low-lying eigenmodes. 


III. CHIRAL CONDENSATE IN QUARK-ANTIQUARK SYSTEM 

In the presence of color charges, there appears a flux-tube of chromo-electric helds, which 
has been observed on the lattice by measuring the spatial distribution of the held strength 
tensor [SHS|- In this section, we investigate the spatial distribution of the local chiral con¬ 
densate qq{x) (|^ around the static color sources. Previously, a related analysis has been 
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Wilson loop time 



FIG. 3: Schematic picture of the flux-tube measurement. Static quark and antiquark are located 
at {R/2,0) and {—R/2,0) on the XF-plane. 

made, but on a single color source, i.e. a Polyakov line [22H21], or at finite temperature 
where the flux-tube is expected to be suppressed [25] . 


A. Partial restoration of the chiral symmetry in the flux-tube 


We investigate the spatial distribution of the local chiral condensate qq{x) around the 
static color sources by calculating a correlation 


_ {qq(S)W(R,T)) ^ 

{qq[x))w = --(99), 


( 10 ) 


{W(R,T)) 

where W{R, T) denotes a Wilson-loop of size R x T. It represents a pair of static quark 
and anti-quark separated by a distance R. The origin of the coordinate is chosen at the 
center of the loop, which stretches along the A-axis, and the Y and Z-axes correspond to 
the transverse directions. Figure shows a schematic picture of the measurement. 

As mentioned in the previous section, we truncate the sum over the eigenmodes in Q at 
the N-th. eigenvalue and denote the corresponding local condensate as qq^^\x). In the hnal 


7 








analysis, we chose N = 160, after confirming that the resnlt is nnchanged once a snfficient 
nnmber of low-lying modes are inclnded. 

Reflecting the nltraviolet divergences of the scalar density operator, the expectation valne 
of (x) contains qnadratic and logarithmic divergences. The strong qnadratic divergence 
is associated with a mixing with the identity operator and has the form mq/a?. Becanse of 
the exact chiral symmetry of the overlap-Dirac operator, the strongest divergence of 1/a^ is 
absent and the leading term is of 1/a? and proportional to Since the truncation at a 
hxed mode number N can be considered as a certain regularization scheme, the regularized 
operator qq^^\x) can be parametrized as 


qqW = ^ cf^rriq/a^ + 


m„ 


( 11 ) 


with the operator for which the power divergences are subtracted. The second and 

third terms represent a mixing with the identity operator; the mass dimension is compen¬ 
sated by irtq/a? and m^, respectively. These coefficients and can be obtained by 
htting the vacuum expectation value {qq^^'^) as a function of the valence quark mass irig 


j. When the correlation with the Wilson-loop is considered as in (10), the contribution 
from the identity operator with the divergent coefficient c^^'^niq/a? -|- cancels on the 

right hand side and the measurement is free from the power divergences. 

Figure l^shows the spatial distribution of {qq^^\x))w on the XR-plane with a separation 
R = 8. The location of color sources are shown by the circles. In order to improve the signal 
of the Wilson loop, we apply the APE smearing for the spatial link-variables, and the 
temporal extent is hxed at T = 4 for which the ground state becomes dominant. In this 
plot, qq{x) is set at t = 0, and the valence quark mass is nig = 0.015a“^. 

In order to improve the signal, the lattice data are averaged over space-time. Namely, 
assuming the translational invariance of the expectation value, we shift the whole system 
including the Wilson-loop and the local chiral condensate and take an average. This can 
be done without additional computational cost to solve quark propagators by using the 
low-lying eigenmodes. This is one of the advantages of the construction (|^. 

As Figure [^demonstrates, there appears a tube-like structure between the color sources, 
where the change of the condensate becomes positive, i.e. {qq^^\x))w > 0. It means that 
the magnitude of the chiral condensate is reduced between the color charges, since (qq) is 
negative in the vacuum. 




FIG. 4: Spatial profile of the local chiral condensate ix))w around a Wilson loop W{R,T) 

with R = 8 and T = 4. The position of color sources are at {X, Y) = (4, 0) and (—4,0), which are 
shown in the plot by white circles. 


Peak structures at the position of the charges are shown in the flux-tube measurements 
due to the strong enhancement of the action/energy density around the color charges. 
In terms of the low-mode truncated local chiral condensate shown in Fig. |^no such charac¬ 
teristic structures around the color charges can be observed. The absence of peaks will be 
discussed later. 

The remaining logarithmic divergence in can be canceled by taking a ratio 


r(x] = — 






( 12 ) 


where is obtained by htting the vacuum expectation value (qq) to (11) as a function 

of the valence quark mass niq. As there are no remaining ultraviolet divergences, the ratio 
r{x) has a proper continuum limit. Hereafter, we mainly use this quantity to quantitatively 
estimate the restoration of chiral symmetry. 

Figure]^ shows the ratio r{x) for the separation between the color sources hxed at i? = 8. 
The plots Fig. (a) and Fig. (b) correspond to the cross-sections of Figure along the 
A-axis and the transverse H-axis. The location of color sources is shown by black dots in 

Fig- 13(a) • 

These plots provide a quantitative measure of the reduction of the chiral condensate. 
The region where the chiral condensate is reduced forms a structure that resembles the 
color flux-tube. In other words, chiral symmetry is partially restored inside the flux-tube. 
The restoration becomes stronger around the center of the flux, which is about 20% when 


R = 
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(a) cross section at F = 0 


(b) cross section at X = 0 
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FIG. 5: Ratio of the chiral condensate r(x) for a separation R = 8. The color sources are separated 
along the X axis and set at {X,Y) = (4,0) and (—4,0). The plots show the cross section (a) along 
the X-axis at T = 0 and (b) along T-axis at X = 0. 


1.0 

0.9 

0.8 

0.7 


• 

• m' 

n m 

m 

a 

m 

m 

m 

m 

m 

m 

m 

n m 

m 

R = 

8, T=0 -a— 









* 

* 

* 


* 

R 

= 8, X=0 


-6 -4 -2 0 2 4 6 

Y 


The close relationship between the flnx-tnbe is snggested in Fig. We compare the 
cross-section of both chiral condensate and the action density defined by ^ 

with a cutoff on the mode number {p^^'>{x))w around the Wilson loop using low-lying 160 
eigenmodes. The latter is calculated inserting the action density p{x) in place of qq{x) 


in Eq. (10), which is used for the flux-tube measurement [3H5]. In order to compare the 
prohle, both quantities are normalized to unity at the origin. Apart from their normalization 
coefficients, the spatial profile of the chiral condensate shows a good agreement with UV 
Dirac mode truncated action density. As mentioned above, the action density is strongly 
enhanced around the color charges as reported in However, both UV hltered densities 

do not have such structures. Our conclusion is that such peak mainly comes from ultraviolet 
divergent part and thus cannot be seen in Fig. (a) within our cutoff scale. 

Figure shows the same plot as Fig. but with different values of N, the number of 
eigenmodes included in the sum (|^. As expected from the construction that cancels the ul¬ 
traviolet divergences, there is no signihcant difference between N = 120 and 240. Our choice 
N = 160 is therefore sufficiently conservative to estimate the local chiral condensate inside 
the tube. Up to the largest eigenmode in our calculation at X = 240, we have conhrmed 
such saturation for other quantities considered in this paper except for the magnitude of the 
action density p^^^(T') and topological charge density q[^J{x). The value of these quantities 
strongly depends on the cutoff scale Amax as expected from the definition in Eqs. (fsl) and 
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(a) cross section at F = 0 


(b) cross section at X = 0 




X Y 


FIG. 6; The spatial profile of both local chiral condensate and UV Dirac mode 

truncated action density {p^^\x))w around the color sources with a separation R = 8 using 160 
low-lying eigenmodes. For a comparison of their shape, both quantities are normalized at the 
origin. 


(a) cross section at F = 0 





(b) cross section at X = 0 



Y 


FIG. 7: Same as Fig. but with different numbers of eigenmodes included: N = 120, 160, 200, 
and 240. 


(|^. However, the spatial profile of both {p^^\x))w and are rather stable, and 

there are no signatnre of peaks within onr trnncation as in Fig. (a). 

The partial restoration of chiral symmetry is in accordance with the chiral bag model 
pictnre for the qnark-antiqnark system m In the naive bag-model, chiral symmetry is 
completely restored inside the bag, while Fig. [^snggests a smooth bonndary with a rednced 
bnt non-zero condensate inside the bag. 
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FIG. 8: Chiral condensate ratio r(x) along the X-axis. The results with increasing separation R: 
i? = 4, 8 and 10. Color sources are located at (i?/2,0) and {—R/2,0). 

B. Chiral symmetry restoration as a function of the separation 

Next we study the chiral symmetry restoration depending on the separation of the color 
sources. 

Figure 1^ compares the cross section of r(x) along the X-axis with i? = 4, 8 and 10. By 
increasing the separation, we observe that the region of partial restoration stretches between 
the color sources, which are located at X = R/2 and —R/2. This supports the picture of 
the tube structure. 

The magnitude of the reduction increases with R. For instance, at the origin, the re¬ 
duction of about 15% at i? = 4 grows up to about 25% at i? = 10. Beyond R = 10, the 
statistical signal becomes much worse, and the effect of spatial boundary would become 
important as R approaches L/2. 

In Figure we plot the value of the ratio at the center r(0), where the magnitude 
becomes minimum, as a function of R. As the separation R increases, the ratio of the chiral 
condensate decreases monotonically until the maximum distance we could explore. At larger 
distances, the effect of string breaking should manifest itself in dynamical QCD, and the 
local chiral condensate would stop decreasing. As far as we can observe, the reduction of 
chiral condensate inside the color flux-tube is of the size of 20-25% at the distance of 1 fm, 
assuming that the string breaking does not occur in this scale [28], since it is difficult to 
observe the breaking state using the Wilson loop as a color source. 

By increasing the separation between color sources, the thickness of the flux is expected 
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FIG. 9: Ratio at the center of the flux r(0) as a function of the separation R. 




FIG. 10: Gross section of the ratio along the R-axis with different separations R and cut points 

X. 


to grow logarithmically as a function of its length [29l |30] . Such behavior has indeed been 
observed in quenched lattice QCD calculations [31] (see also [32] for a study at finite temper¬ 
ature), for which one can increase the statistics more easily. In this work we try to observe 
the thickness through the chiral condensate ratio r{x). 


Figure 10 shows the cross-section of r{x) along the F-axis for some combinations of R 


and X. It is clear that for a fixed X the flux is thicker when the separation R is larger. More 
interestingly, the curve for i? = 8 at X = 0 almost coincides with that for i? = 10 at X = 2. 
Similarly, the curve for i? = 4 at X = 0 coincides with that for i? = 6 at X = 2. Due to 
the reflection symmetry, these behaviors are also observed at X = —2. This indicates that 
the thickness of the flux is highly correlated with the magnitude of the reduction. We also 
note that these corresponding cross-sections have the same distance from the color charge. 
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R X 

r(0) 

f fi a 

X^/dof 

10 

0 

0.757(10) 

1.54(7) 0.66(11) 2.3(0.9) 

0.40 

10 

1 

0.762(10) 

1.46(7) 0.72(14) 2.8(1.2) 

0.46 

10 

2 

0.778(10) 

1.26(6) 0.71(14) 2.5(1.1) 

0.66 

10 

3 

0.805(9) 

1.00(6) 0.75(18) 2.5(1.3) 

1.00 

8 

0 

0.786(5) 

1.11(3) 0.71(7) 2.2(0.5) 

0.43 

8 

1 

0.792(5) 

1.08(4) 0.72(8) 2.3(0.6) 

0.25 

8 

2 

0.813(5) 

0.93(3) 0.75(11) 2.5(0.8) 

0.58 

8 

3 

0.855(5) 

0.69(3) 0.83(17) 3.0(1.3) 

1.42 

6 

0 

0.815(3) 

0.89(2) 0.66(4) 1.7(2) 

0.38 

6 

1 

0.827(3) 

0.81(2) 0.65(4) 1.6(2) 

1.01 

6 

2 

0.865(3) 

0.65(2) 0.61(4) 1.4(2) 

2.36 


TABLE I: Fit results for r{Y) to (13) for each R and X, together with per degrees of freedom. 


Note that /i and a are in lattice units. The condensate ratio at the center r(0) is listed as well. 


which is given by R/2 - |X| for |X| < R/2. 

In fact, snch coincidence is expected from an effective string model [33l |3l]. According 
to that model, the ratio r{Y) is written as 

( 13 ) 

a A'i(q!) 

where Kq{x) and Ki{x) are modihed Bessel functions. The parameter /i has a physical 
interpretation as the inverse penetration length of the flux from the perpendicular direction, 
and a is the thickness of the core. The parameter r represents the strength of the condensate 
reduction. 


The function (13) reproduces the lattice data quite well, as shown by the curves in Fig. 10 


The £t results are summarized in Table The penetration length is in the range of 1/p ~ 
1.0-1.6, which corresponds to 0.11-0.18 fm in physical units. The core size a is 1.4-3.6 in 
lattice units, and is in the range 0.15-0.4 fm. We observe an increase of a as i? increases 
while X is hxed at zero, but with the large statistical error we are not able to claim the 
clear evidence of the string fattening. 
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IV. CHIRAL CONDENSATE IN THE THREE-QUARK SYSTEM 


A. Partial Restoration of chiral symmetry in the 3Q-system 


Next we consider a system consisting of three color charges that represents a baryon 
system, that we call the 3Q-system in this paper. 

Using the path-ordered product Uk = Hrfe along a path Tk, the 3Q Wilson-loop is 

given by 

W"3Q = IsabcSa'b'c'Ur'UfUf, (14) 

which is made color-singlet by the totally anti-symmetric tensor Sabc of color indices a, b, 
and c [35l |36]. Similar to the QQ-system, the spatial distribution of the chiral condensate 
for the 3Q-system is measured as 




( 16 ) 


with the 3Q Wilson loop IUsq. The ratio of the chiral condensate in the 3Q-system r 3 Q(x), 
for which the ultraviolet divergences cancel, is then constructed by 


r3Q{x) = 


(16) 


(gg(™bt))(pp3Q) • 

Figure [TT] shows a schematic picture of the construction of the 3Q Wilson loop IF 3 Q from 
the Wilson lines Uk- For simplicity, we use an isosceles right triangle conhguration of the 


color charges on the W-plane, and the coordinate is set as in Fig. In this case, the 
junction point of the three flux tubes (the Fermat point) corresponds to the origin [351136] . 
The measurement of the local chiral condensate qq{x) is done at a hxed time slice. The 
low-mode truncation number N, the temporal extension T and other measurement setups 
are the same as in the QQ-system. 

shows the ratio tsq^x) with color sources at {X,Y) = (6,0), (0,6) and (0,0) 


Figure 
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denoted by circles in the plot. As shown in Fig. the magnitude of the chiral condensate is 
reduced among the color sources, which indicates the partial restoration of chiral symmetry 
inside the 3Q-system. Similar to the QQ-system in Sec. Ill, there appear no peaks at the 
color charges within our truncation scale. We note that the characteristic Y -type flux is not 
clearly seen in this plot, probably because the thickness of flux is comparable to the color 
source separation. Because of the statistical noise, we are not able to repeat the calculations 
increasing the quark separations. 
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3Q Wilson Loop 



FIG. 11: A schematic picture of the construction for three-quark system with an isosceles right 
triangle configuration. Three Wilson lines Uk correspond to the static color sources. 





FIG. 12: Condensate ratio r 3 Q{x) with the color sources at Qi = (6,0), Q 2 = (0,6), and 
Q 3 = (0,0) on the AT-plane. 


Like in the QQ-system, the magnitude of the restoration depends on the separation of 
shows the cross-section of the ratio r 3 Q(x) along the line of X = Y 


the sources. Figure 
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with the color sources at {X,Y) = {R,0), {0,R) and (0,0). In this setup, the measurement 
goes through one color charge and the center of mass of the system. By comparing the data 
for i? = 3 and 6, we hnd that the reduction is more substantial for R = 6, which is similar 
to the QQ-system (see Fig. [^. The reduction of the local chiral condensate becomes larger 
with the size of the loop, and take its minimum value at around the center of gravity. With 
R = 6, the reduction is about 30%, which is also similar to that of the QQ-system. 
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FIG. 13: Chiral condensate ratio r 3 Q(af) along the line of X = Y with the color sources at 
{R, 0), (0, R) and (0,0) on the XY-plane with R = 3 and 6. 

B. Partial restoration at finite density 


Finally, using the observed modification of the local chiral condensate around the color 
sources, we estimate the size of the partial restoration of chiral symmetry in finite density 
QCD. We consider the system of fixed number of baryons in a finite volume box, so that 
the baryon number density p is N^/L^, where Nj, is the number of baryons and is the 
spatial volume. As a toy example we take A'";, = 1 and replace the baryon by the 3Q 
Wilson-loop. This only gives a crude approximation of the realistic system, but given the 
difficulty of simulating QCD at finite chemical potential it may provide a useful clue to the 
understanding of the finite density QCD. 

The net change of the condensate under such system is estimated by the spatial average 
of the condensate ratio r^q^x): 


{m)p 

(99)0 


1 

1 ? 


L3 


E-3Q(f), 


(17) 


where {qq)p is the condensate at the finite baryon number density p = l/L^. We use 
two lattice volumes, = 16^ and 24^, which correspond to (16a)“3 ~ 0.18 fm“^ and 
(24a)“3 ~ 0.05 fm“^, respectively. The lO^ lattice roughly corresponds to the normal nuclear 
density po — 0.18 fm“^. 


Figure 14 shows {qq)p/{qq)o as a function of 1/L^. The two symbols correspond to the 
different configurations of the color sources, i.e. (0,0), (i?, 0) and (0,7?) with R = 3 and 
6 on the XD-plane. The solid lines are the results of a linear fit with fixed value of 1 at 
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FIG. 14: Reduction of the chiral condensate at finite density measured by {qq)p/{qq)o- The 
estimates for a given configuration of color sources, i.e. at (0,0) and (0, R) on the XY- 
plane, with R = 3 and 6 are shown. 


l/L^ = 0. The linear dependence from unity at p = 0 simply means that there is a hnite 
region where the chiral condensate is reduced from its vacuum value. Since the region gets 
larger with increasing R, the slope for larger R is steeper. 

In our setup, the reduction of the chiral condensate at the normal nuclear density is only 
~ 5%, which is much smaller than the phenomenological model estimate of the order of 30% 
133. Our estimate, however, assumes a hxed spatial size of the baryon which is smaller than 
the realistic nucleon. For instance, the mean root square radius of our setup in Fig. 14 is 


0.44 fm when R = 6, while the charge radius of proton is 0.88 fm. As the restoration of 
chiral condensate is stronger for larger separation, this suggests that {qq)p in realistic hnite 
density QCD could be substantially lower than our estimate. 


V. SUMMARY 

The Dirac eigenmodes carry the full information of the background gauge held. Indeed, 
having the complete set of the eigenvalue and eigenvectors, one can reconstruct the held- 
strength tensor F^y{x) at any point x. They therefore oher an interesting way of hltering out 
the ultraviolet modes and investigating the low-energy dynamics of QCD by only using the 
low-lying eigenmodes upon reconstruction. This is a sound and well-dehned regularization 
method of quantum held theory. 

We use this method to investigate the spatial prohle of the chiral condensate under 
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the presence of external sonrces. On the lattices generated with 2+1 flavors of dynamical 
overlap fermions, we calculate the low-lying eigenvalues and associated eigenvectors of the 
overlap-Dirac operator, and use them to reconstruct the chiral condensate locally. Then, it 
is straightforward to measure its correlation with the external color sources set up to model 
the QQ and 3Q systems. 

We hud that the local chiral condensate shows a structure interpreted as a color flux-tube 
between the QQ color sources, in which the condensate decreases signihcantly. It indicates 
a partial restoration of chiral symmetry inside the flux-tube and suggests that it happens 
also inside hadrons. The spatial prohle is consistent with a string model of the conhnement 
potential, giving another support for the presence of the color flux-tube. 

We perform a similar measurement in the 3Q system, which is new as far as we have 
noticed. It again shows the partial restoration of chiral symmetry among the color sources. 
The reduction of condensate is about 30% for the separation between the color sources of ~ 
1 fm. It can be used to estimate the chiral condensate in the hnite density system. 

The method developed in this work may easily be applied for the study of hnite tempera¬ 
ture QCD, where Polyakov loops can be used for a static color source. Since the eigenmodes 
can be applied to dehne various charges, such as the axial charge density, the quark number 
density, and the topological charge, it may provide an interesting alternative to measure 
their spatial distribution. 

Acknowledgements 

The lattice QCD calculations have been done on SR16000 at High Energy Accelerator 
Research Organization (KEK) under a support of its Large Scale Simulation Program (No. 
13/14-05). This work is supported in part by the Grant-in-Aid of the Japanese Ministry of 
Education (No. 25287046 and 26247043), and the SPIRE (Strategic Program for Innovative 
REsearch) Field 5 project. 


[1] T. Banks and A. Casher, Nucl. Phys. B 169, 103 (1980). 

[2] K. Fukushima and T. Hatsuda, Kept. Prog. Phys. 74, 014001 (2011) |arXiv:1005.4814 [hep- 
ph]]. 


19 



[3] G. S. Bali, K. Schilling and C. Schlichter, Phys. Rev. D 51, 5165 (1995) |hep-lat/9409005|. 

[4] R. W. Haymaker, V. Singh, Y. C. Peng and J. Wosiek, Phys. Rev. D 53, 389 (1996) |hep- 
lat/9406021,. 

[5] P. Cea and L. Cosmai, Phys. Rev. D 52, 5152 (1995) |hep-lat/9504008,. 

[6] C. Gattringer, Phys. Rev. Lett. 88, 221601 (2002) |hep-lat/0202002i. 

[7] S. Aoki, T. W. Ghiu, G. Gossu, X. Feng, H. Fukaya, S. Hashimoto, T. H. Hsieh and T. Kaneko 
et al, PTEP 2012, 01A106 (2012). 

[8] H. Neuberger, Phys. Lett. B 417, 141 (1998) |hep-lat/9707022|. 

[9] H. Neuberger, Phys. Lett. B 427, 353 (1998) |hep-lat/9801031|. 

[10] H. Fukaya et al. [JLQCD Collaboration], Phys. Rev. Lett. 98, 172001 (2007) jhep-lat/0702003 . 

[11] H. Fukaya et al. [TWQCD Collaboration], Phys. Rev. D 76, 054503 (2007) |arXiv:0705.3322 
[hep-lat]]. 

[12] H. Fukaya et al. [JLQCD Collaboration], Phys. Rev. Lett. 104, 122002 (2010) [Erratum-ibid. 
105, 159901 (2010)] [arXiv:0911.5555 [hep-lat]]. 

[13] H. Fukaya et al. [JLQCD and TWQCD Collaborations], Phys. Rev. D 83, 074501 (2011) 
arXiv: 1012.4052 [hep-lat]]. 

[14] T. Iritani, G. Cossu and S. Hashimoto, PoS LATTICE 2013, 376 (2014) |arXiv:1311.0218 
[hep-lat]]. 

[15] T. Iritani, G. Cossu and S. Hashimoto, PoS Hadron 2013, 159 (2013) jarXiv: 1401.4293 [hep- 
lat]]. 

[16] T. Iritani, G. Cossu and S. Hashimoto, PoS (Lattice 2014) 338 |arXiv: 1412.2322 [hep-lat]]. 

[17] F. Niedermayer, Nucl. Phys. Proc. Suppl. 73, 105 (1999) |hep-lat/98100261. 

[18] H. Fukaya et al. [JLQCD Collaboration], Phys. Rev. D 74, 094505 (2006) |hep-lat/0607020|. 

[19] S. Aoki, H. Fukaya, S. Hashimoto and T. Onogi, Phys. Rev. D 76, 054508 (2007) 
arXiv:0707.0396 [hep-lat]]. 

[20] E.-M. Ilgenfritz, K. Roller, Y. Roma, G. Schierholz, T. Streuer and V. Weinberg, Phys. Rev. 
D 76, 034506 (2007) |arXiv:0705.0018 [hep-lat]]. 

[21] E.-M. Ilgenfritz, D. Leinweber, P. Moran, R. Roller, G. Schierholz and V. Weinberg, Phys. 
Rev. D 77, 074502 (2008) [Erratum-ibid. D 77, 099902 (2008)] [arXiv:080I.1725 [hep-lat]]. 

[22] W. Eeilmair, M. Faber and H. Markum, Phys. Rev. D 39, 1409 (1989). 

[23] W. Sakuler, W. Burger, M. Faber, H. Markum, M. Muller, P. De Forcrand, A. Nakamura and 


20 


I. O. Stamatescu, Phys. Lett. B 276, 155 (1992). 

[24] M. Faber, M. Schaler and H. Gausterer, Phys. Lett. B 317, 409 (1993). 

[25] S. Chagdaa and E. Laermann, PoS LAT 2007, 172 (2007). 

[26] J. Noaki, T. W. Chiu, H. Fukaya, S. Hashimoto, H. Matsufuru, T. Onogi, E. Shintani and 
N. Yamada, Phys. Rev. D 81, 034502 (2010) |arXiv:0907.2751 [hep-lat]]. 

[27] A. Hosaka and H. Toki, Phys. Kept. 277, 65 (1996). 

[28] G. S. Bali et al. [SESAM Collaboration], Phys. Rev. D 71, 114513 (2005) |hep-lat/0505012|. 

[29] A. Hasenfratz, E. Hasenfratz and P. Hasenfratz, Nucl. Phys. B 180 (1981) 353. 

[30] M. Liischer, G. Munster and P. Weisz, Nucl. Phys. B 180 (1981) 1. 

[31] N. Cardoso, M. Cardoso and P. Bicudo, Phys. Rev. D 88, 054504 (2013) |arXiv: 1302.3633 
[hep-lat]]. 

[32] A. S. Bakry, D. B. Leinweber and A. G. Williams, Phys. Rev. D 85, 034504 (2012) 
arXiv: 1011.1380 [hep-lat]]. 

[33] P. Cea, L. Cosmai and A. Papa, Phys. Rev. D 86, 054501 (2012) jarXiv: 1208.1362 [hep-lat]]. 

[34] D. E. Kharzeev and F. Loshaj, Phys. Rev. D 90, 037501 (2014) |arXiv: 1404.7746 [hep-ph]]. 

[35] T. T. Takahashi, H. Matsufuru, Y. Nemoto and H. Suganuma, Phys. Rev. Lett. 86, 18 (2001) 
[hep-lat/00060051. 

[36] T. T. Takahashi, H. Suganuma, Y. Nemoto and H. Matsufuru, Phys. Rev. D 65, 114509 (2002) 
[hep-lat/02040111. 

[37] R. S. Hayano and T. Hatsuda, Rev. Mod. Phys. 82, 2949 (2010) |arXiv:0812.1702 [nucl-ex]]. 


21 


